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^ : Abstract 






We use Monte Carlo simulations to study the static and dynamical properties 
O ' of a Potts glass with infinite range Gaussian distributed exchange interactions 

for a broad range of temperature and system size up to A^ = 2560 spins. The 
results are compatible with a critical divergence of the relaxation time r at the 
theoretically predicted dynamical transition temperature To, r oc {T — Tu)^ 
with A ~ 2. For finite N a further power law at T = To is found, t(T = 
Td) oc N^ with z* « 1.5 and for T > To dynamical finite-size scaling seems 
to hold. The order parameter distribution P{q) is qualitatively compatible 
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'^ . with the scenario of a first order glass transition as predicted from one-step 

^ [ replica symmetry breaking schemes. 
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f^^ • Developing a theory of the glass transition of a fluid from its Hamiltonian within first 

^ . principles statistical mechanics is still a formidable problem [0^- While some researchers 

attribute glassy freezing to the (hypothetical) vanishing of the configurational entropy |^ 



O 
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O ! at the "Kauzmann temperature" Tx (which is lower than the experimental |I| glass 

^ ' transition temperature Tg), others emphasize the dynamical transition at the critical tem- 

c^ ■ perature Tc of mode coupling-theory [Q from the ergodic fluid to a non ergodic state. Since, 

for atomic systems, Tc > Tg, this frozen phase can have only a finite lifetime and is thought 
to decay by thermally activated (so-called "hopping") processes. 
C I Recently evidence has been given p,^ that these two seemingly different scenarios could 

§ I both result as two complementary aspects of the same unifying theory 0]. In view of the 

questions that still exist on the various theoretical approaches, it is valuable to have exactly 
solvable models that exhibit a similar behavior: a dynamical transition at a temperature To 
and a static first order glass transition at a temperature Tq < T^- One of these models is the 



5^ ! p-state infinite range Potts glass with p > 4 [§-[T3|, where at Tq a static (Edwards- Anderson 



type |T^|T^) spin glass order parameter qea appears discontinuously. However at Tq there 
is neither a latent heat nor a divergence of the static spin glass susceptibility Xsg- The 
latter would diverge only at an extrapolated spinodal temperature Tg < Tq, see Fig. 1. The 
dynamical behavior of the spin autocorrelation function C{t) for T > To is described by 
the same type of equations [0,0 as found in mode-coupling theory ||^. Thus this model 



seems indeed to have many properties in common with structural glasses. Apart from being 
a possible prototype model for the structural glass transition, the Potts glass can also be 



considered as a simplified model for an anisotropic orientational glass: e.g. a six-state Potts 
glass may be a reasonable description of a diluted cubic molecular crystal where diatomic 



molecules can align only along the p = 6 face diagonals [|1^ . 

In order to understand this model in more detail, we have performed extensive Monte 
Carlo simulations. Of course these simulations were done for systems with finite size and 
thus all the transitions in Fig. 1 are rounded. While in the past finite size effects at normal 
first and second order transitions have been studied extensively |T^,|TB[, very little is known 
about finite size effects at dynamical transitions. Thus it is hoped that the present work will 
be useful for the proper analysis of simulations of realistic models for the structural glass 



transition as well |T9|. In addition, finite mean-field systems may also have some similarities 
with systems having finite interaction range [^: the relaxation time r is large but finite 
for T < Td as well, unlike the behavior in the thermodynamic limit N ^ oo (shown in Fig. 
1), since energy barriers between "basins" are finite for T < To pO[. Thus the mean field 
Potts glass for finite N should be a good model for elucidating the physics of glassy systems 
in general. 

We now give some technical details about our simulations. We study the Hamiltonian 



1 ^ 



(1) 



where each spin cr^ G {l,...,p} interacts with all the other spins. The interactions Jij 
are taken from a Gaussian distribution with mean Jq = (3 — p)/{N — 1) and variance 
AJ = (N — 1)~^/^. We choose p = 10 (note that for p = 2 the standard Ising spin glass 
results [0]). For this choice of A J it can be shown P that T^ = 1 |^. We simulate the 
system sizes iV=160, 320, 640, 1280 and 2560, and the number of independent samples is 
500 for A^=160 and T > 1, else 100, with the exception of iV=2560 where only 50 samples 
where used. The equilibrium dynamics is studied by means of the standard Metropolis 
algorithm |1l5|-[T7[]. In order to produce equilibrated configurations and to study the static 
properties we used again the Metropolis algorithm at high temperature (T ^ T^), while for 
lower temperatures we used the parallel tempering method p2| , ^ . 



Fig. 2 shows typical data for the autocorrelation function C{t) of the Potts spins ||2^ 
at two temperatures and various system sizes. While at very high temperature, such as 
T = 1.8, finite size effects are completely negligible, at T = To they are unexpectedly 
pronounced. Whereas in the thermodynamic limit one expects at the dynamical transition 
C(t -^ oo) = q^A, for A^ ^ 1280 hardly any indication of the development of the plateau (as 
expected from Fig. 1) is seen. This behavior is strikingly different from the behavior found 
in analogous autocorrelation functions for atomistic models of structural glasses, where a 
pronounced plateau occurs already at temperatures a few percent above Tc (= Tq in our 
model) [0 and finite size effects are hardly detectable for a comparable range of A^. While 
for N ^ oo the dynamics of the present model for T ^ T^ is described by mode-coupling 
theory ||12|,|13|, the assertion that infinite range models for finite A^ resemble the behavior of 



real structural glasses |2^ seems, in view of Fig. 2, doubtful to us, at least for the system 
sizes presently accessible. 

In order to quantify the slowing down of the dynamics as a function of the temperature 
and system size we define a relaxation time r by C{t) = 0.2. This particular choice was 
made since we must have C{t) < qEA(T = To). Our results indicate that, for T = Tq, r 



behaves like a power law t oc N 
analytical results 



finite size scaling hypothesis 
of N and T for T > T^ , 



with z* ^ 1.5 (see inset in Fig. 3). It is known from the 

131 that, in the limit N —^ oo, t{T) oc (T/Td — 1)~^. The dynamical 

T7|,p!8|,p5| assumes a generalized behavior for r as a function 



r = N'>{N{T/Td - 1)^/'*} for A^ ^ oo and {T/Td - 1) 







(2) 



with the scaling function f(^ -^ oo) oc ^^^* to recover the proper thermodynamic limit. As 
can be seen from our data in Fig. 3, this ansatz is satisfied by the Potts glass in the vicinity 
of the dynamical transition. Our choice of A/z* Ri 1.3 in Fig. 3 implies A Ri 2, which is 
similar to values found for atomistic models of the glass transition [|19| and compatible with 
a direct extrapolation A^ — ^ cxo of the data for t{T, N) for T > Td |p6 |. 

We mention that Eq. (El) has a well-based theoretical foundation for second order phase 



transitions 



^,^], where in Fig. 1 Tq = T^ = Tg would coincide. However, for a dynami- 
cal transition which has no associated diverging static susceptibility this relation is purely 
phenomenological. For second order mean-field spin glass transitions one has Eq. (^ with 
A/z* = 2/3mf + IMF = 3, where Pmf = 1,7mf = 1 are the static mean-field exponents of 
the spin glass order parameter and the susceptibility, respectively [^. Since A = 2 ||T5[ 
one has z* = 2/3, consistent with expectation |2^ and simulations ^U|. Obviously, in view 
of the systematic deviations from scaling that are still visible in Fig. 3, our estimates for A 
and z* clearly are somewhat tentative only. 

We conclude by discussing some results about some static properties of this model. As 
usual in numerical simulations of spin glasses, the order parameter is a measure of the overlap 
between microscopic states visited by two different replicas (same realization of disorder, but 
different dynamics) of the system. For Ising spin glasses, this is just the number of spins in 
the same state divided by the total number of spins. In the Potts glass, as well as in vector 
and quadrupolar spin glasses, the overlap between spins belonging to different rephcas is a 
tensorial quantity. Therefore we define an overlap order parameter which is invariant under 
global rotations of the spins 



15,16|: 
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p-1 . N 

Y.ir'^Y with q'" = j^Y^^t^^h. 



(3) 
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where a and /? are the replica indices and S'f q, are the components of the spins in the simplex 
representation |]2^. The static spin glass susceptibility is defined as xsg = ^[{llip)]- Our 
Monte Carlo results (Fig. 4) are indeed compatible with the expected behavior for static 
quantities: the spin glass susceptibility remains finite at T = T^, and the order parameter 
distribution P{q) develops a two peak structure not by splitting off a single peak when the 
temperature is lowered, as would be common for second-order transitions [|l5H l8|. Instead 
the second peak grows near q = qea continuously gaining weight for low temperatures. This 
temperature dependence is consistent with a one-step replica symmetry-breaking scenario 







In summary we have shown that the finite size rounding of the dynamical transition in 
the p = 10 infinite range Potts glass is compatible with a finite size scaling hypothesis, and 
that the relaxation time at the dynamical transition scales like a power law, r oc A^^*, with 
z* ~ 1.5. The static spin glass susceptibility converges to a finite results at To, as expected. 



In contrast to atomistic models for the glass transitions, this model allows well equilibrated 
simulations at T = Tp and also much lower temperatures (Fig. 4), at least for A^ in the range 
of a few hundreds. While atomistic models allow to study questions such as the dynamical 
heterogeneity only at T > Tjj (which would not be a relevant temperature region in this 
context for experiments), the present model allows to determine the relaxation behavior of 
various quantities over a much wider temperature range. Hence from the present model we 
expect that a rather complete picture of the glass transition can be obtained, which thus 
should give stimulating insights into the proper analysis of simulations for more realistic 
models too, and help in clarifying the fundamental questions that still remain about the 
nature of glassy freezing in structural glasses. 
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FIG. 1. Qualitative sketch of the mean- field predictions for the p— state Potts glass model 
with p > 4. The spin glass order parameter in thermal equilibrium is positive only for T < Tq 
and jumps to zero discontinuously at T = Tq, where the spin glass susceptibility XSG is finite (for 
T > Tq, xsg follows a Curie- Weiss type relation with an apparent divergence at Tg < Tq). The 
relaxation time diverges already at the dynamical transition temperature Tq. This divergence is 
due to the occurrence of a long lived plateau in the time-dependent spin autocorrelation function 
C{t). 
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FIG. 2. Spin-spin autocorrelation function C{t) for T = 1.8 and for T = To = 1.142 |11], for 
several values of N. The solid line is the theoretical value of the Edward- Anderson order parameter 
1ea{Td) ioi N —> qo |11|. The dashed line locates the value we use to define the relaxation time 
r. 
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FIG. 3. Log-log plot of the scaled relaxation time t/N^* vs. the scaled distance in temperature 
N{T/T£) — 1)^/^* from the dynamical transition temperature, choosing z* = 1.5, A/z* ~ 1.3. The 
inset is a log-log plot of t{T = To) vs. A^. 
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FIG. 4. Inverse of the spin glass susceptibility Xsg plotted vs. T ^. Different symbols corre- 
spond to different values of A^. The number of different samples used is 300 for N = 160, 200 for 
N = 320 and 100 for N = 640 and 1280. Arrows locate Ts,To,Td Q. The inset shows the order 
parameter distribution P{q) vs. q for A^ = 320 and various temperatures. 



